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Abstract — In this paper, a novel technique is proposed to 
address the joint sampling timing acquisition for baseband and 
broadband power-line communication (BB-PLC) systems using 
Orthogonal-Frequency-Division-Multiplexing (OFDM), including 
the sampling phase offset (SPO) and the sampling clock offset 
(SCO). Under pairwise correlation and joint Gaussian assump- 
tion of received signals in frequency domain, an approximated 
form of the log-likelihood function is derived. Instead of a 
high complexity two-dimension grid-search on the likelihood 
function, a five-step method is employed for accurate estimations. 
Several variants are presented in the same framework with 
different complexities. Unlike conventional pilot-assisted schemes 
using the extra phase rotations within one OFDM block, the 
proposed technique turns to the phase rotations between adjacent 
OFDM blocks. Analytical expressions of the variances and biases 
are derived. Extensive simulation results indicate significant 
performance improvements over conventional schemes. Addition- 
ally, effects of several noise models including non-Gaussianity, 
cyclo-stationarity, and temporal correlation are analyzed and 
simulated. Robustness of the proposed technique against violation 
of the joint Gaussian assumption is also verified by simulations. 

Index Terms — OFDM, power-line communication, baseband 
system, sampling phase offset, sampling clock offset, non- 
Gaussian noise 



I. Introduction 

Power-line communication (PLC) is a potential solution 
for long haul transmission, last mile access, and in-building 
connection over low to high voltage power networks. By 
turning virtually each wired device into a target of future 
value-added services, PLC acts as a technological enabler 
spanning internet access, residential or business premises, 
smart grid, and other municipal applications [1 1. Moreover, it 
enjoys the advantages of pervasive power cable infrastructures 
and low implementation cost without rewiring [2|. 

To provide a higher data rate than earlier ultra narrow- 
band PLC (UNB-PLC) and narrowband PLC (NB-PLC) solu- 
tions, broadband PLC (BB-PLC) empowered by multi-carrier 
schemes increases the data rate to 200 Megabits per second 
(Mbps) and the frequency band to 1.8 — 30 MHz (3j. One step 
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further, the new ITU and IEEE standards extend the frequency 
band to 100 MHz and increase the data rate to 500 Mbps 0. 

Originally designed for electricity delivery, power line is 
a harsh and noisy medium, particularly for low power and 
high frequency data transmission [1|. Statistical modeling of 
PLC channels is a technical challenge for communication 
theory approach [0]. In 0, Zimmermann et al. suggested 
an analytical model for the channel transfer function (CTF) 
of PLC, characterized by a small number of parameters for 
frequency below 30 MHz. To broaden its coverage, Tonello 
et al. in [6| proposed a random broadband channel generator 
based on measurement campaign results on nine classes of 
channels in the 2 ~ 100 MHz range. J7), O discussed 
the time-variation of PLC channels. However, for broadband 
transmission, it could be largely neglected, since its period 
(reciprocal to the mains frequency) is much longer than the 
block duration. 

Different types of noise are observable in power lines. A 
previous work [2| suggested the Nakagami-m distribution to 
fit the envelope probability density function (PDF) of time 
domain background noise. Its accuracy was proved in [9| by 
empirical measurements in a 10-kV medium-voltage power 
network. Relevant works can be found in ifTUl . ifTTI . The Class- 
A impulsive noise model was proposed in lfl2l . It has been 
extensively used ever since due to its canonical and closed- 
form PDF Q~3), as well as excellent agreement with measure- 
ments from both natural and man-made noise environments 
fl4l . The power spectral density (PSD) of impulsive noise 
can reach a value of more than 50 dB above background noise 
lfT31 . In lfl6l . Katayama et al. revealed the cyclo-stationarity of 
the noise, synchronous to the mains voltage frequency. Also, 
noise samples are temporally correlated, as shown in ifPTl . 

Orthogonal-Frequency-Division-Multiplexing (OFDM) is 
one option for multi-carrier schemes in BB-PLC. For baseband 
OFDM systems, signals are transmitted without up-down 
conversion onto carrier frequency. The complex conjugate 
property, known as the Hermitian Symmetry Property (HSP) 
1 18 1, must be satisfied in frequency domain after modulation 

im. 

Using Fast Fourier Transform (FFT), OFDM significantly 
enhances the system performance in a dispersive propagation 
environment lEOl . Nevertheless, it is vulnerable to synchro- 
nization imperfections. In this work, two related issues caused 
by a non-ideal sampling clock at the receiver are considered, 
including: 
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1) Sampling phase offset (SPO): It is caused by the mis- 
aligned initial sampling instants between the transmitter 
and the receiver, resulting in an extra phase rotation 
proportional to the tone index within one OFDM block 
(inter-block phase rotations) (181 . It is equivalent to a 
fractional symbol timing offset (STO) l2D . Although 
it could be incorporated into CTF and eliminated by a 
frequency equalizer (FEQ) |fT8l , the residual part should 
be compensated due to oscillator instabilities (22l . 

2) Sampling clock offset (SCO): It is caused by a non- 
synchronized sampling clock between the transmitter 
and the receiver, triggering not only a inter-block phase 
rotation, but also a phase growing linearly for successive 
OFDM blocks (intra-block phase rotations) (TP . 

Both SPO and SCO could be tracked recursively using pilots 
by a delay-locked-loop (DLL) (22l . (23l . The loop operation 
demands an accurate initial timing acquisition. To the authors' 
best knowledge, joint sampling timing estimators are seldom 
discussed in literature. Thus, for comparison, pilot-assisted 
SPO and SCO estimators are investigated: 

• SPO estimators: [24] proposed a minimum mean square 
error (MMSE) estimator for narrowband signals. A low 
complexity estimator was proposed in [25] taking advan- 
tage of the phase rotation between two adjacent pilots. 
(26) presented several improved variants to better utilize 
the phase and power characteristics within one or several 
OFDM blocks. The scheme in (27) further reduced the 
pilot distortion using the amplitude information and thus 
enhanced the performance in a multipath channel. Chan- 
nel state information (CSI) is explicitly required in [24], 
while implicitly in [25|-[27|. 

• SCO estimators: [28] proposed an estimator using the 
phase rotation between the upper and the lower frequency 
band without CSI knowledge. In (29), Morelli et al. 
proposed a maximum likelihood decoupled estimator 
(MLDE) and a reduced-complexity estimator (RCE) with 
superior performance in a multipath channel. CSI could 
be estimated by the maximum likelihood criterion in [29 1. 
The scheme in (30) was initially devised for MIMO- 
OFDM, which was modified to suit a single-antenna 
system in (29). (311 proposed a weighted joint least 
square estimator (W-JLSE) outperforming conventional 
LSE estimators ED. Both ED and ED require CSI. 

In the present work, a novel joint acquisition technique for 
SPO and SCO is proposed. An approximated form of the log- 
likelihood function is derived under pairwise correlation and 
joint Gaussian assumption of received signals in frequency 
domain. To avoid the high complexity two-dimension (2- 
D) grid-search which cannot be decomposed into a one- 
dimension (1-D) one due to mutual dependency, a five-step 
method is proposed. First of all, ancillary phases are obtained 
block-wise, either by a non-data-aided (NDA) estimator or 
by its data-aided (DA) counterpart. Then, phase unwrapping 
restores the actual phases from the wrapped ancillary ones. 
The resultant phases are further embedded into a linear model. 
With such linearity, LSE is used to separate SPO and SCO. 
Bias erasure removes the inherent bias of estimation, followed 



by another LSE which leads to refined estimations. The 
proposed technique highlights the intra-block rotations with 
performance improvements over conventional pilot-assisted 
schemes. Effects of various types of noise and violation of 
the joint Gaussian assumption are analyzed. 

The remainder of this paper is structured as follows: Section 
iHl introduces the signal model. Section [Til] illustrates the 
approximated log-likelihood function, several variants of the 
proposed technique, and analytical expressions for variances 
and biases. Section [IV] discusses the applicability to time 
domain non-Gaussianity, cyclo-stationarity, and temporal cor- 
relation. Section [V]presents simulation results. Finally, Section 
[Vl] concludes the paper. 

II. Signal Model 

For Q OFDM blocks in frequency domain with N sub- 
carriers in each block, the output signal s(t) at the transmitter 
on baseband is |23| 



q=o fee/c 



(1) 



where 



/,\ 1 i2irk( 



Ng + qNs 



' u(t - qN s T sam ) (2) 



is the subcarrier pulse, JC the locations of the data subcarriers, 
X q} k the equi-probable data modulated on the fc-th subcarrier 
in the q-th OFDM block and X q> k € T; T denotes the finite 
alphabet of constellation. N s = N + N g is the length of an 
OFDM block, N g the length of the cyclic prefix (CP), and 
T S am the sampling interval at the transmitter. The rectangular 
function u(t - qN s T sam ) is 



u(t) 

The multipath channel is 

h{t,r) = 



1 0<t<N s T sam , 
otherwise 



L-1 

E 

1=0 



hi(t)6(r- n ) 



(3) 



(4) 



where L denotes the total number of taps, {/i;(i)}i=o,i,... .l-i 
the real channel gains, {t;}; = o ! i.... .l-i the time delay of each 
path, and <£(•) the delta function. Transmitting s(t) through the 
multipath channel (@) yields 



r(t) = s(t) * h(t, t) + w(t) = V hi(t)s(t - n) + w{t) (5) 



L-1 

E 

1=0 



'The following notations are used: (■)* as the complex conjugation; (-) T 
the transposition; (-) H the Hermitian transposition; E{-} the expectation; 
Var {■} the variance; 5R{-} the real part of its argument; the imaginary 

part of its argument; Z the set of all integers; Iff the N X N the identity 
matrix; J jv the N X N anti-diagonal matrix; Ojvj x N 2 me Ni x N% all zero 
matrix; CA r (n(fJ.), c(<r)) the complex Gaussian distribution with mean fi 
(or mean vector fi) and variance a (or covariance matrix <r); Tr{-} the trace 
of a square matrix; max[-] the maximum value of the arguments; diag{X} 
the diagonal matrix with entries given by vector X; Ex me autocorrelation 
matrix of X defined as E {XX fl }; the pseudocorrelation matrix of X 
defined as E{XX T }; the (i, j)-th entry of a matrix; Other notations 

would be either self-evident or explained therein. 
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where * stands for linear convolution, and w(t) is the noise 
with an unknown distribution and correlation. 

At the receiver, sampling r(t) at time instants nT^ am gives 
L-i 

r ( nT Lm) = h l( nT Lm) S ( nT Lm - T l)+ W ( nT 'sam) ( 6 ) 

1=0 

where T' sam is the sampling interval at the receiver. Af- 
ter removing CP, the N samples in the q-th OFDM block 
is denoted as r q = [r g>0 , r q>1 , ■ ■ ■ ,r g>n ,-- - ,r q> N-i] where 
r q ,n = r ((n + qN s + N g )T' sam ). However, a timing offset 
£ exists due to timing imperfections, which consists of the 
integer part Int{£} (caused by the misalignment between 
the correct starting position of FFT window and estimated 
position) and the fractional part Frac{£} (caused by the 
SPO). Without loss of generality, a left-shift of FFT window 
is assumed with Int{£} G \—{L — 1),0]. Thus, r q>n — 
r((n + qN s +N g +Z)Ti am ). 

The channel is assumed to be quasi-static within one OFDM 
block's duration. Therefore, hi(q) denotes the real channel 
gain of the Z-th path in the q-th OFDM block. After demodu- 
lation of the q-th received OFDM block r q by FFT , the k-th 
subcarrier is given by ll20l 

A 

JV-1 

x q , k ,H q , k ,e^ k ' k '^n(k'(i + n) - k) +w q , k 



where R q = [R qfi R qA 



Ra 



k'=Q 
k'^k 



(7) 



where 



JV-l 



W 9lk = $>((« + Q N > + N 9 + OZam) 

(8) 

Hi)*) = M-TT (9) 



<f>(q,k,k',£,r)) = 



N sin 

n(N - l)(fc'(l +v)~k) + 2tt(1 + 7y)£fc' 



N 

2nqN s r]k' + 2TcN g r)k' 
N 



(10) 
(11) 



Here, H q>k = ^f^ 1 hi{q)e^ fZ7T ^rt^ is the CTF of channel 



-j2n 



1, and n = {T' sam - T sam )/T sam the SCO. The diagram 
illustrating SPO and SCO is drawn in Fig. Q] The extra phase 
rotation Q q /. £ „ can be approximated by 

QftM.u" -^-+7T??fc+ — (12) 

for large N and small rj. H(r]k) is the amplitude attenuation 
which approaches 1 for small 77 [32 1 . lCI^ yn - q . k is the inter- 
carrier interference (ICI) caused by £ and rj. 

For convenience, lCl^ tV - qyk is neglected; its effect is demon- 
strated in Section [V] Thus, (01 can be written into a compact 
matrix form below: 



q,N-l] T , *<? 

[0,N - 1]}, H„ 



diag{& 7 ' e «' fc - s - , », k 
diag{H g }, H«j = [H q , H q> \ ■ ■ ■ H q ,N-i] J ', X g = 
[X 9i0 • • • X q>N _x] T , W q 4 [VF g ,o W 9) i • • • W^v-i] T 
Here, we make two assumptions: 

(Al) n q X q ~CAf(0 Nxl ,<T 2 x <T 2 H I N ). 

(A2) W 9 ~ CA/"(0jvxi, Oypliv) an d independent from 

HqXq. 

Therefore, R 9 is a y'o/nf Gaussian vector distributed as 
CA/"(0jvxi, (<4ff| + cr^)Ijv). On the other hand, the psew- 
docorrelation matrix (PCM) of H 9 X 9 and W 9 can be written 
into 



iV 



[O, otherwise 
'a^y, i+j = N 
0, otherwise 



(14) 



(15) 



using the proprieties and HSPs of H g X 9 and W g respec- 
tively. For instance, the propriety of W q>k indicates that 
E {W qtk W q , k >} = 0ifk^N-k', while the HSP implies that 

W q,k = W q N-k and so E{W q ,kW q , k >} = Ow if and only if 
k = N — k'\ 
Remarks: 

In practice, neither X q k nor H q k is Gaussian distributed: 
X q}k is taken from a finite alphabet with equal probability 
and thus is uniformly distributed, while the statistical model 
in J6l is adopted for H q:k . The closed-form distribution of 
H f/ X g is hard to pursue in practice. Even when both X q , k and 
H q:k are Gaussian distributed, their product H q:k X q:k follows 
the Gaussian product distribution which is non-Gaussian. 
(Al) and (A2) significantly simplify our analysis by making 
Rq,k ~ CM{Q,<j\afj + Robustness of the proposed 

schemes against non-Gaussianity of R 9 is demonstrated via 
simulations in Section fVl 



■W„ 



III. Proposed estimators 

A. Approximated Log-likelihood Function 

An OFDM system with N n null subcarriers at both ends 
of the spectrum as the guard bands is considered here. The 
universal set containing all the N subcarriers is denoted as 
Xjj, which can be decomposed into three subsets with index 
given by 

. X = {k\k = 0, k = N/2, fceZ}=l„ + U Xq 
. Zi = {k\k = {N -N n )/2 + v,0 < v < (N n /2)-l,v€ 
Z} U {k\k = (N/2) +v,l<v<(N n + l)/2 + 1, v £ 
Z} = 1+ U Xf 
. 1 2 = {k\l < k < (N — N n )/2 — l,k £ Z}U{fc|(N- 
N n )/2 + 2< k< 7Y-l,JfeeZ}=2j UX^ 
For BB-PLC baseband system, Xq delivers random real-valued 
data, which are often nullified since direct current (DC) and 
Nyquist frequency subcarriers are generally discarded. X\ con- 
tains the null subcarriers, while subcarriers in I2 transmit data. 
Each subset can be divided into the left and the right half. The 
cardinalities for Xq ~ X2 are: Card{Zo} = 2, Card{2i} = 
(13) N n , Card{I 2 } = A^ u 4 N - N n - 2 
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Fig. 1. Diagram for SPO, SCO, Int{£}, and Frac{£} for the q-th OFDM block. 
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Subcarrier Index 



Fig. 2. Decomposition into Xq = Xq U X , X\ = X^~ U I 1 , and I2 = X^ U I 2 



- 1 fcta» + 2 jv - 2 a' — 1 



For the g-th block, an alternative sequence can be formu- 
lated by taking the conjugation of the right half of R g : 



R 



q,0 



R. 



q,N/2 



-1 R 



q,N/2 



R 



q,N-l . 



(16) 



N/2 



N/2 



Now, the autocorrelation matrix of could be written into 



U W 
(JV-l)xl 







lx(iV-l) 
B 



(17) 



where the sub-matrix B is shown in ( fTSI l on the next page. 
Ctv„+i represents the (N n + 1) x (JV n + 1) bi-diagonal matrix 
with unity diagonal and anti-diagonal entries^, and 9 .e = 
&q,k,C,r, + Qq,N-k,i,r, = 2tt£ + T) (irN + 2irN g + 2TTqN s ') Vfc 
where 9 = [£ r]) T . B n is closely tied with the noise 
structure since B„ = [Sw]i :JV -i,i : iv-i + [ s w]i:JV-i,i:JV-i 
where [■}n 1 :N 2 ,n 3 :N 4 represents the sub-matrix spanning rows 
Ni ~ N2 and columns A3 ~ N4. Discussions of B„ under 
various scenarios are given in Section ITVl 

Arranging observations of Q OFDM blocks into a vector 
R' = [Rq R' x Rj • • • Rq_-J, after manipulations reported 
in Appendix [Aj the approximated log-likelihood function con- 
ditioned on 9 under (Al) and (A2) is 

Q-i 

A(R'|0) = Const. + 2 ^ ^{Mq,^-^*- 9 } (19) 
9=0 keX + 

X 1 {q,k)=R q<k R qtN - k (20) 

Const, is irrelevant to estimation and can be neglected. 
A(R'|0) coherently accumulates ^R.{Xi(q, fe)e _J s- 9 } across 
all Q OFDM blocks and subcarriers in I 2 . 



-For instance, 



B. Joint Estimation of £ and r/ 

Joint estimation of 9 using JT9l takes the form 

9 = argmaxA(R'|0) 
e 

,Q-i \ 
^argmaxl ^ ^ R{Ai(g, fc)e^' '- e } ) 



(21) 



where 9 is the estimation of 9. (|2~TT > requires a 2-D grid search 
on all possible 9 for the global maxima. The complexity 
involved can be formidably high and cannot be pursued in 
practice ||29l . One seemingly plausible solution is to decouple 
the original 2-D problem into a 1-D one as in [29 ] . However, 
it is not possible since estimations of £ and r\ are mutually 
dependent: estimating £ by fixing 77 introduces bias into £ and 
vice versa. A five-step method is employed to overcome this 
issue: 

1) Ancillary Phases Estimation: Obtaining ancillary 
phases Q qt o,Vq e [0, Q — 1] by either a NDA or a DA 
estimator. For the q-th OFDM block, 



Q q ,e = 



ar g{ E fce i+ A i(<7> fc )j> 
ar S^ Efc e i+ h(q, k)X 2 {q, k) 



NDA 
DA 



(22) 

where X 2 (q, k) is introduced to highlight the modulation 
effect and frequency selectivity as 

X 2 (q,k) = [2 + SNR- 1 (q,fc)] _1 (23) 

where SNR(<j, k) = igs^jga*!! . This is explained in 
more details in Appendix lAl For notational convenience, 
A denotes either NDA or DA. The DA suits the closed- 
loop system with feedback to compute X 2 {q, k). On the 
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(J X CJ H L N u /2 

(A r „ + l)xA' u /2 



B 







JV„/2x(JV„ + l) 



<7 X f7 H) 



'N u /2 















JV u /2 Ojv u /2x(7V„ + 1) 



(iV„+l)xiV u /2 



a W^-N u /2 

(jV„ + l)xAT„/2 
a W^N u /2 











A r „/2x(JV„ + l) 
N u /2x(N„+l) 







(AT„+l)xiV„/2 



(18) 



other hand, the NDA, with a similar form to that in 
reduces the complexity without calculating \2(q, k). 

2) Phase Unwrapping: Since Q q +i,e - @ q ,o = Const., 
{Qq.e, Q = 0, 1, • • ■ , Q — 1} lie on the same line. How- 
ever, arg in (|22| > folds into [— ir, +ir] and breaks the 
line into several discrete segments, known as the phase 
ambiguity. To deal with this issue, phase unwrapping is 
used to detect the discontinuities in phases sequentially 
and restore the actual phases. More specifically, once the 
criterion 

l©£T* - e&l > ^ (24) 

is satisfied, ±2pir is added to the phases starting from 
e , until (1241 is met again (and ±2(p + 1)tt is added 
afterwards). <p e [0, 2] is the tolerance level of phase 
unwrapping, p denotes the number of times that (|24| | is 
met before 0^ e . The sign in (|24T > is determined by the 
direction of phase alternation. 

3) Linearization and Least Square Estimation: The un- 
wrapped phase is denoted as 0^ e . They could be com- 
bined into the vector 0^. Now, assuming correct phase 
unwrapping across all the Q OFDM blocks, can be 
embedded into a linear model shown as 

@ A =m A +V (25) 

The Q x 2 observation matrix E can be expressed by 

-I T 



E 



2tt 
Dq 



2?r 
Di 



2- 



2n 
Dq-i 



(26) 



where D q = ttN 



2nNg + 2irqN s and V is the self- 
noise with mean vector /xf; SNR and covariance matrix 

^0,SNR 



, expressed by 



0.SNR 



-,0,SNR 



sin(6g,e 



cos(6 



SNR 



< q < Q - 1 



(27) 



1 - cos(e g . e ) 
diag ^ -TV x SNR 



< q < Q - 1 



(28) 



which are derived under (A2) and (A3) , (A4) . (A3) , 
(A4) are shown below: 
(A3) fiat fading (L = 1) with real channel gain 
ho(q) with zero mean and variance a\. Thus, 
in frequency domain, EjiJ^fe} = and 
E{\H q , k \ 2 }=a 2 H . 
(A4) Constant modulus modulated signal satisfying 
E^fel^OandEjlA^fcl 2 }^. 



Notice that, (A3) and (A4) differ from (Al) : they 
merely impose the finite first moment and second moment 
on X q .k and H q ^ without the Gaussian assumption on 
HqX 9 . 3 = jfi is the ratio between usable pairwise 
correlations to the useful FFT duration, and SNR = 
yUi^jL where v = jf- compensates the loss of discarding 
unused subcarriers. In lack of 6, it is impossible to pre- 
compute /x v and S-y . An ordinary least square 
(OLS) estimation yields the initial A : 



Oa = (E T E)" 1 E T 0^ = 6 A + V 



(29) 



where V = (E T E) _1 E T V is the 2 x 1 error vector 
associated with V by a linear transformation (E T E) _1 E. 
4) Bias Erasure: With A and assuming perfect a priori 
knowledge of SNR, both /x^ NR and £^ SNR can be 

replaced by /x v and 2j v . Bias erasure (BE) 
gives 



f 0,4, SNR • 

of fi^' in 







0, 



0.A.SNR 



E9 A + V 



(30) 



where M ^' SNR =^' SNR -Mv' SNR andS^' SNR = 

5) Re-estimation: Updating the estimation using the 
weighted least square estimation after bias erasure (WLS- 
BE) by 



WLS-BE 



= (E 2 



V<0A, SNR 
2J V , 



E)" 1 



E 1 



-,6» A .SNR 



0. 



(31) 



or OLS after bias erasure (OLS-BE) ignoring the weight- 
ing factors 







A 



OLS-BE 



(E t E) _1 E t 0.a 



(32) 



The performance loss of (l32t comparing with (BTt could 



be neglected, since X 



0.A.SNR 
V 



approaches a diagonal 



matrix with nearly identical entries as cos(0 9i e) varies 
slowly block-wise given small -q. 



C. Variance and Bias of Estimation 

Var {£}, Var {r}}, Bias{£}, and Bias{^} before BE with 
perfect phase unwrapping can be computed straightfor- 
wardly from (f25} under (A2) , (A3), and (A4) . The de- 
tails are shown in Appendix [B] Analytical expressions 
of Var {£},Var {7]},Bias{^}, and Bias{f/} are given by 
(l33Tl.d34Tl.d5Fb.j36t respectively. 
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Var{£} = 



Var {rj} ■■ 
Bias{£} 



2g 2 (2Q + 1)(Q + 1) + 2<?(Q + 1)(4Q - 1) + (2Q + 1)(2Q - 1) 



E?=o 1 cos(0,, e ) 



4tt 2 7V(1 + g) 2 Q{Q 2 - 1)SSNR 
(Q-l)[l + 4g + 4Q(l + 5 )] -6g[ fl + Q(l + 5 )] 



4tt 2 (1 + <?) 2 iVQ 2 (Q + 1) 2 (Q - 1) 2 SSNR 
3 S?=o 1 9cos(e 9 , fl )(2g-Q + 1) 2 



7r 2 7V 3 (l + g) 2 Q{Q 2 - 1)SSNR tt 2 (1 + ,g) 2 AT 3 Q 2 (Q + X )2(g _ i)2 SSNR 



sin(6g,, 



O-i 

" S C0S ( e «. e ) + SNR 

O-i 



(Q - 1)[3 + 6 5 + 2(1 + ,g)(2Q - 1)] - 6g[l + 2g + (1 + ,g)(Q - 1)] 



27r(l + fl )Q(Q + l)(Q-l) 



Bias{?7} = ^ 



sin(e g: e) 



^ cos(e g , e ) + SNR 



3(Q - 1) - 6g 



7r(l + g)iVQ(Q + l)(Q-l) 



(33) 
(34) 
(35) 

(36) 



Here, g = -jf- is the ratio between CP and useful FFT 
duration with details in Appendix [B] The mean square errors 
(MSEs) are 

MSE {£} = Var {£} + [Bias{£}] 2 (37) 
MSE {fj} = Var {rj} + [Bias{^}] 2 (38) 

a SNR 

With perfect BE, /x v , = 0<gxi> an d consequently the bi- 
ases in ( f35T > and (|36] l vanish which makes MSE{£} = Var {£} 
and MSE{^} = Var{rj}. 
Remarks: 

(i): The first terms of (l33~t and (l34l are written as Var{£}i 
and Var{?7}i. It is easy to prove that max Var{£} = 

2Var{£}i and max Var {rj} = 2 Var {rj} i. They provide a 
good criterion on the choice of Q and N. For instance, when 
N is fixed, to warrant the performance of MSE {^{Target and 
MSE {fj]Target under a certain SNR assuming perfect BE, Q 
should be chosen according to ( f39l > shown on the next page, 
for Q > 1 by setting MSE {^Target = max Var {£} and 

MSE {rj} T , 



arget 



iax Var {rj} ; ft denotes the SNR gap 
considering the possible loss of accuracy of (l33l and (|34V 

(ii): Obviously, A2(g, k) in (l23l reduces to SNR (g, fc) in low 
SNR, and becomes 0.5 in high SNR. Thus, the advantage of 
DA is more significant in low SNR . 

(Hi): For correct phase unwrapping in absence of noise, two 
conditions must be satisfied: 

(CI) -7T < 27r£ + nr)N + 2-KN g r\ < +tt for correct £ 
estimation 

(C2) — 2ir < 2irN s ri < 2tt for correct r\ estimation 

Note that, only Frac{f} is estimable, since Int{£} results in 
an imperceptible 27rlnt{£} rotation in 6 g .e; the only impact 
of Int{£} on Frac{£} is an increased ICI^ )Tj;9) fc. Thus, £ is 
treated equivalently as Frac{£} and bounded in [—0.5,0.5]. 
The hexagon denoting the solution area is plotted in Fig. [3] 
with area S = ?( N + N J ' w hich shrinks in presence of noise. 

(iv): In practice, \ 2 (q : k) (for the DA scheme) and SNR 
(for both DA and NDA schemes) need to be calculated. Firstly, 



A2(g, k) in (1231 is computed as 

h(<hk) = [2 + SNR-!(g,fc) 
_ ( \x^\STkl 



(40) 



X qt k is either 



where SNR-^fc) 

X q ,k (for pilots or null subcarriers), or the output of decision 
device X q k- H q: k stands for the estimated CTF With ifTSl. 
(Tyy and SNR can be estimated by 



r2 — 



Ai(g,fc) 



SNR 



X U H 



■w 




where 



9 y^Q-i 



\R. 



q,k\ 



'X 



'X U H 



'W 



QNn 



(41) 
(42) 

(43) 



Now, SNR is inserted back into 



and 



/^v . Replacing in d30l > with /x v ^ 



to yield 
R leads 

to &a =&a Mv A SNR ■ Usin S OLS-BE in (02) gives the 
practical estimation 0^| pg , termed as the practical scheme 
(PS) given by: 



6a\ = (E T E)- 1 E T & 



IPS 



(44) 



IV. Applicability to Different Types of Noise 

A. Class-A Impulsive Noise 

The PDF for time domain real Class-A noise random 
variable (r.v.) w is shown in Table U [12|. T denotes the 
FFT operation. T = <7g/af is the power ratio between 
the Gaussian component with variance a 2 and the additional 
man-made noise component with variance a 2 . The Class- 
A model combines both components and its PDF can be 
regarded as a weighted sum of infinitely many Gaussian PDFs 
with an increasing variance a 2 . Impulsive index A is the 
product of the average rate of impulsive noise and the mean 
duration of a typical impulse. The noise is more impulsive 
for a small A, whereas for A 00, the PDF becomes 
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Q 



^ (MSE{e}Tar get ^ 2 S(SNR + l ,^(l[SE[/^r,„, / ./7r(I- ; /)-E(SXR + <>] 



(39) 



V 



~A7 



m = 



N+N g 



N+N s 



6 = -0.5 



6 



m = 1-2g 

N+2N„ 



'' 2 N+2N q 



0.5 



Fig. 3. The region of £ and r; for perfect phase unwrapping in a noiseless environment. The coordinates for Ai ~ Ag are: Ai : (—0.5, (N + N g ) 1 ), A2 

{-N g (2(N + N g ))' 1 ,(N + N g )- 1 ), A 3 : (0.5,0), A 4 : (0.5, (TV + N,)- 1 ), A5 : (N g (2(N + 2V S )) _1 , — (JV + -/V 9 ) _1 ), A 6 : (-0.5,0). 

TABLE I 

Statistical information for Class-A and Nakagami-m distribution 



Noise 



Parameters 



PDF (time domain) 



PDF (frequency domain) 



Class-A 



A, T 



/W = £r=0 7§j-exp(-^) 



W = T{w} 



2 V/A+T 

a T 



Nakagami-m 



m, n 



/(«) 



2m"'(n) z 



r(m)n™ 
n 2 



n = e{« 2 } 

= \/2mcos(C), C ~ U{-tt, +tt) 



= Jf{5R{«;}} 

w ~ ca^(o, n) 



a Gaussian distribution 1141 . After FFT, TE in frequency 
domain approaches Gaussianity, particularly for a large A, 
which is accurate for p > 3 fl3l . 

B. Nakagami-m Background Noise 

The envelope PDF of time domain Nakagami-m background 
noise r.v. u is given in Table U [10], where m is the fading 
figure controlling the severity of the amplitude fading [34 1 
and m > 0.5. T( ) is the Gamma function, and f2 the second 
moment of u [35 1 . For the corresponding complex random 
noise w with an envelope following the Nakagami-m PDF 
(|tu| = u), the axis-PDF for its real part = ucos(£) 

is shown in power series in [2| and in closed-form in [ 1 1 ; 
£ is uniformly distributed in [— n, +7r]. Yet, the closed-form 
distribution for the corresponding frequency domain r.v. W is 
hard to pursue. As a rule of thumb, W can be modeled as a 
Gaussian r.v. with zero mean and variance ft, by virtue of the 
central limit theorem and the uniformly distributed ( [2|. An 
additional y/2 is inserted to compensate the loss of ignoring 
3{ui}. 

Remarks: 

(i): In a strict sense, m differs in different frequency ranges. 



For high frequencies, m ~ 1, while m < 1 for low frequencies 
[2J. For simplicity, a common m for all subcarriers is assumed. 

(ii): A good tool to evaluate the goodness-of-fit is the % 2 - 
test. In this paper, the D'Agostino's x 2 -test is adopted with 
a significance level a — 0.05, which is based on the trans- 
formation of the third and fourth-order statistics (skewness 
and kurtosis) on data samples 11361 . T is set as {1,0.1,0.01} 
with respect to light, moderate, and heavy impulsiveness, while 
A <G {0.01,0.1,1}. For the Nakagami-m background noise, 
m £ {0.5,1.0,2.0,10}. The kurtosis and skewness for the 
Gaussian case are 3 and respectively. 

Independent and identically distributed (i.i.d) noise samples 
are generated using the corresponding PDFs. The averaged 
in-phase p-value (probability that the distribution cannot be 
rejected as Gaussian) on usable subcarriers with N = 512 
is shown in Table HjH. Also, the kurtosis and skewness are 
tabulated. Similar results can be obtained for the quadrature. 
For most cases, near-Gaussian performances can be achieved. 
Thus, Sw, S^y, and B„ in ( fT8b are invariant. The Gaussianity 

3 For the Nakagami-m background noise, the DC subcarrier could reach a 
considerably high value which would destroy the Gaussianity. Fortunately, it 
is usually discarded in most OFDM systems. 
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TABLE II 

KURTOSIS, SKEWNESS, AND P-VALUE OF THE D'AGOSTINO'S X 2 -TEST FOR CLASS-A AND NAKAGAMI-M NOISE {N ■ 



512, a = 0.05) 



Noise 


Parameters 


Kurt . ( i n - ph as e ) 


Skew.(in-phnse) 


p-value(in-phase) 






A = 0.01 


9 ftl 


Q 10 \y in - 4 

y . i. a x J.U 


n KQOQ 




1 = U.01 


A = 0.1 


9 CM 


7 91 v i n — 4 


u.yooz 






A = 1 


9 Q8 

z.yo 


o in v 1 — 3 


u.yozi 


Class-A 




A — n m 


9 RQ 


9 7 v i n ~~ 3 

— £t. I X J.U 




T = 0.1 


A = 0.1 


2.94 


— 1.3 x 10~ 4 


0.9544 






A = 1 


2.98 


8.75 x 1CT 4 


0.9498 






A = 0.01 


2.90 


6.78 x 10~ 4 


0.9435 




T = 1 


A = 0.1 


2.98 


1.90 x 10~ 3 


0.9502 






A = 1 


2.98 


-8.1 x 10" 4 


0.9461 




m = 


= 0.5 


2.98 


-8.3 x 10- 4 


0.9154 


Nakagami-m 


m — 


- 1.0 


2.99 


-1.4 x 10~ 3 


0.9173 


m = 


: 2.0 


2.98 


-2.4 x 10" 4 


0.9116 




m = 


10.0 


2.98 


-1.1 x 10" 3 


0.9112 



degrades for A = 0.01 with T = 0.01 or T = 0.1. 

C. Effect of Temporal Correlation and Cyclo-stationarity 

For colored noise, the zero vectors in B„ are replaced by 
non-zero values, and both the diagonal and the anti-diagonal 
elements of B„ are not uniform and this degrades the esti- 
mation. The cyclo-stationarity could be analyzed in a similar 
way. The uniformity of the (anti-)diagonal elements in B„ 
is preserved which entails no performance loss. Meanwhile, 
non-zero elements appear in the original zero vectors in B n . 
As the zero vectors in B„ are not involved in the estimation, 
they would not degrade the performance. 

V. Simulation 

An OFDM system with the following parameters is consid- 
ered, unless otherwise mentioned: 

1) Sampling interval: T sam = 10ns (sampling frequency 
100 MHz) 

2) Number of subcarriers: N = 512 

3) Guard interval: N g = 64 

4) OFDM blocks: Q = 10 

5) Null subcarriers: N n = 64 

6) Modulation: 16-PSK (M = 16, a\ = 1) 

7) Channel: The real channel taps are generated from the 
model proposed in J6). Channel classes 1 ~ 9 with 
different channel capacities are selected with respect to 
the probability profile {0.0349, 0.1678, 0.1818, 0.1188, 
0.1188, 0.1258, 0.0979, 0.0769, 0.0769}. The number of 
paths L follows a Poisson distribution. Also, the channel 
length is tailored to satisfy L < = 32. Other 
parameters are chosen according to |]6]. 

8) Impulsive parameters: A,T G {0.01,0.1,1} x 
{0.01,0.1,1} 

9) Nakagami-m fading figure: m G {0.5, 1.0, 2.0, 10.0} 

10) SPO: £ = 0.1 

11) SCO: r?= 1 x 10~ 5 

12) Phase unwrapping tolerance level: (p = 1 in (1241 

We generate i.i.d impulsive Class-A noise samples using 
the toolbox in |37|, while i.i.d Nakagami-m background noise 
samples by the square root of Gamma distributed r.v.. MSE{£} 
and MSE{^} are chosen as the criteria given by MSE{£} = 



10 



10 



10 



10 



10 



10 




15 20 25 

SNR(dB) 

Fig. 4. Comparison of MSE{r;} with PS and other schemes. 



E {|£-£| 2 } and MSE{^} = E{|^-^| 2 }. The achievable 
variances after BE are plotted as Var {£} and Var {rj} using 
the first terms of (l3~3~i and (l34l as benchmarks. 

As far as the colored noise, the 1/ f" model is used with 
(3 6 {0.72, 0.337} for generation of correlated noise samples 
[38 1. The condition Tr{B„} = (N— l)a^y is imposed to keep 
the total noise power invariant. 

For the cyclo-stationary noise, similar to the noise model 
proposed in ifTBI . the instantaneous variance is set as 



E{ W (t)}4 a p (t)=A a (*) B m a (^) 



(45) 



where Tag iS the main voltage frequency (60 Hz here). A q 
is chosen according to (|46T > to make the average energy in 
the g-th block [(no + (<? - l)N)T sam , (no + qN - l)T sam ] 
invariant: 



t=(n + (q-l)N)T 3 , 



Na 2 w 



,{n +qN-l)T st 
1/2 



E 



no+qN-l 
n=n +(q-l)N 



jin2 



(46) 
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-A-DA.16-PSK, PS 
— 0— NDA, 16-PSK, PS 
—A- DA, 16-QAM, PS 
NDA, 16-QAM, PS 




Fig. 5. Comparison of MSE{£} with PS and other schemes. 




Fig. 6. Simulation time under different OFDM block number Q. 

A. Comparison of £ and rj estimation performance using PS 
with other schemes 

Fig. g] and Fig. |5] demonstrate MSE{??} and MSE{cf} 
performances of several schemes versus SNR in additive 
white Gaussian noise (AWGN). All available subcarriers are 
treated as pilots (N p = N - N n = 448). Perfect knowledge 
of CSI is assumed for ||24)-||27], (29)-GT] and DA-PS, while 
unknown for NDA-PS and |28). For Oberli's estimator ll30l . 
£ = or otherwise rj would be biased. Likewise, 77 = for 
||24)-||27). RCE is selected from (29), and the Algorithm 2 
out of the four schemes in ||26l . Note that, ||26l and J27| are 
originally designed for estimation of both Int{£} and Frac{£} 
in a multipath channel without knowledge of CSI. However, 
simulations indicate that CSI should be obtained to guarantee 
the accuracy of Frac{f }. 

PS in (02) is simulated under: (SI) DA, 16-PSK; (S2) 
NDA, 16-PSK; (S3) DA, 16-QAM; (S4) NDA, 16-QAM. 
For the performance of rj, (SI) slightly outperforms (S2) 
when SNR < 20 dB. Advantage of DA is not significant for 
moderate SNR since \%{q, k) w Const.. The performance 
gap is more significant in low SNR . For (S3) and (S4) , the 




-y— MSE{?}. DA.OLS 
-y- - MSE{f}, DA.OLS 

MSE{i?},DA,WLS-BE 
MSE{f},DA.WLS-BE 



- MSE{j7}, NDA.PS, K = 0.1 
-< ■ - MSE{f}, NDA.PS, n — 0.1 
— M — MSE{i?}, DA,PS, re — 0.1 
- * ■ - MSE{f }, DA.PS, «= 0.1 





MSE{i7}.NDA,OLS 
MSE{f},NDA,OLS 
i— MSE{?}, NDA.WLS-BE 
■A-- MSE{f},NDA,WLS-BE 
MSE{i)}, NDA.OLS-BE 
■■>-■ MSE{f},NDA,OLS-BE 



15 20 25 
SNR(dB) 

Fig. 7. Different variants of MSE{»7} and MSE{§} estimation including 
OLS, WLS-BE, OLS-BE, and PS with DA and NDA. 



performance gap is enlarged due to the modulation effect of 
QAM . In the ICI-dominant region, (S3) and (S4) outperform 
(SI) and (S2) . 

Similar conclusions can be drawn on £ except that (SI) 
slightly outperforms (S3) when SNR £ [4, 15] dB. According 
to the later analysis, the OLS using DA outperforms NDA 

SNR 

in low to moderate SNR. However, /x^ , derived under 
(A4) , departs from the case of QAM . This is more significant 
for £ according to d35l l and (l36l l. Therefore, an inaccurate OLS 
might in reverse compensate the BE inaccuracy and leads to 
better performance after BE. 

As observed in Fig. Hand Fig. both DA-PS and NDA-PS 
outperform Il24l - ||3TI . It is even more remarkable considering 
that the pilot-assisted schemes exploit all usable subcarriers 
and perfect CSI. 

The averaged simulation times of DA-PS, NDA-PS, and 
[ 24]— [ 31] are plotted in Fig. [6] The complexity of the proposed 
PS schemes are moderate among all the schemes. 



B. Comparison of the proposed DA/NDA variants 

Fig. [7] compares MSE performances of OLS d29), WLS- 
BE (ED, OLS-BE CH, PS (HO) using either DA or NDA. For 
the DA-PS scheme, to simulate imperfect channel estimation 
and therefore imperfect \2(q,k), a N x 1 complex Gaussian 
random vector is introduced which satisfies HSP This vector, 



denoted with J„ 
and variance fc =" 



\Hq.k 



Jq.i ' ' ' Jq,N—i] > nas zero mean 
for each OFDM block, independent 
of H q . The imperfect CTF, denoted as H 9 , is expressed by 



N 



H g = y/l - K 2 llq + KJ q 



(47) 



where k represents the channel estimation accuracy and is 
set at 0.1 here. Before BE, DA-OLS outperforms NDA- 
OLS in low SNR. BE is important for £, while rj does not 
improve after BE mainly due to imperfect phase unwrapping 
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01. A 1 




01. A 1 


-A- r 

■-O-T 


0-1. A 001 
0-1. A 0.1 




0-1. A 1 


—A- T 


l.A 0-01 


— G- .T= 


1,A=0.1 


— *- T= 


1,A=1 




15 20 25 30 35 40 
SNR(dB) 

Fig. 8. MSE{»7} under Class-A and Nakagami-m background noise using 
NDA and PS. 




Fig. 9. MSE{»7} under cyclo-stationarity and temporal correlation using 
NDA and PS. 



not considered in derivations. The performances of WLS-BE, 
OLS-BE and PS are similar. DA outperforms NDA in all 
situations for low SNR . k does not result in any degradation, 
since \2(q, k) is not distorted under (|47| >. 



C. MSE{?y} in Different Noise Models 

Fig. [8] shows MSE{r/} under Class-A impulsiveness and 
Nakagami-m background noise using NDA-PS; in the former 
case, SNR is defined as the background SNR as va\a 2 H / a 2 g 
to stress the effect of T. 

For the Class-A case, a wide margin exists between 
MSE{^} with different T. The performance degradation is 
nearly 7.5 dB when T = 1 => T = 0.1, and 10 dB when 
T = 0.1 =>■ T = 0.01. The gap between AWGN and 
T = 1,A = 1 is 3 dB. A = 0.01 degrades the performances 
for T = 0.01/0.1 due to the reduction of Gaussianity shown in 
Table EH reaching 3.4 dB and 3.3 dB for T — 0.1 and T = 0.01 
respectively; a simple remedy is to increase the FFT size N. 



10" 



10 



<g 10" 

a 

< 10 



ID1: Uniform X, Tonello H 
jD2: (Al) 



— 0— MSE{S;},DA,PS,D1 
--€>-MSE{f},DA,PS,Dl 
— 0— MSEft}, DA,PS,D2 
--0--MSE{f},DA,PS,D2 




— B— MSE{^},NDA,PS,D1 
--B-MSE{f},NDA,PS.Dl 
— #— MSE{tj},NDA,PS,D2 
■-•■-MSE{f},NDA,PS,D2 



15 20 25 
SNR(dB) 

Fig. 10. Compaiison of MSE{5y} and MSE{£} under two statistical 
assumptions. 



The losses for A = 0.1 and A = 0.01, T = 1 stem from the 
increased variance of the numerical result Var {W} = ^-^-a^ 
given a small A. In other words, the numerical variance in the 
Class-A case in Table Ull converges only in the mean sense for 
finite size simulations. 

As far as the Nakagami-m background noise, performances 
under different m are indistinguishable and approach the 
performance under AWGN, since the noise variance is de- 
termined by f2 and irrelevant to m. 

Fig.|9]shows MSE{?7} under cyclo-stationarity and temporal 
correlation using NDA-PS. The cyclo-stationarity entails no 
performance loss, since the uniformity of the (anti-)diagonal 
entries of B„ is preserved as mentioned in Section |IV] On 
the other hand, the temporal correlation degrades the perfor- 
mances, which is more pronounced for /3 = 0.72. 

The same conclusions can be made on DA and MSE{£} 
performances which are omitted here due to space limitation. 

D. MSE{r?} and MSE{£} With Different Statistical Models 
of Signal and Channel 

Fig.[T0]compares MSE-j^} and MSE{£} with PS using DA 
and NDA under 

(Dl) : Uniform distributed X g and H 9 in [6] (Kurtosis: 
8.8120) 

(D2) : Gaussian distributed H g X g under (Al) (Kurtosis: 
2.9741) 

For (D2) , samples of H 9 X 9 are generated directly from 
CAf(0jvxi, Oxcrffljv). Interestingly, MSE {f}} with DA under 
(Dl) slightly outperforms that under (D2) , while (D2) 
slightly outperforms (Dl) for SNR G [2,15] dB in terms 
of MSE{£}. The kurtosis of (Dl) is significantly higher 
than that of (D2) , signifying a sharper peak and longer tail 
distribution. Thus, the probability for a considerably high 
amplitude H^X^ under (Dl) outnumbers that of (D2) , which 
enhances the profit of \2(q,k) using DA. For MSE{£}, as 
stated in lV-AI (D2) gains from the reverse compensation. The 
same performances can be obtained for NDA under (Dl) and 
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(D2) . The small performance gap implies the robustness of 
the proposed technique against violation of (Al) and thus the 
joint Gaussian assumption of R g . 

VI. Concluding remarks 

A novel five-step joint timing acquisition technique is 
proposed to estimate SPO and SCO. An approximation form 
of the log-likelihood function is derived under the pairwise 
correlation and Gaussian assumption on the received data in 
frequency domain. By using the intra-block phase rotations, 
the proposed technique achieves significant performance im- 
provements over conventional pilot-assisted estimators relying 
on inter-block phase rotations. Analytical expressions for the 
variances and biases are derived. Effects of several noise 
models including non-Gaussianity, cyclo-stationarity, temporal 
correlation are discussed. The robustness of the proposed 
technique against violation of the joint Gaussian assumption 
is verified by simulations. 
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Appendix A 

Derivation of the approximated log-likelihood 
function in ([T9li under (Al) and (A2) 

Using ( TToT l under (Al) and (A2) in SectionUll R g is a joint 
Gaussian vector. Similar to ll40l . the log-likelihood function is 
derived as 

,Q-l 

A(R'|0)=log( II II W,k,K,N-k) 

(A.48) 

f(R' q , k ,R' q , N -k) = ^ de t [SG ] ex P (-GPgI^G*) 

(A.49) 



where G 4 [R' q ^R^ N _ k ]. S G = E {G H G} is given by 



°X°H + ° 



W 



'X U H 



■1°w 



'X U H^1 U W U X U H^ U W 

(A.50) 

For convenience, 7 is introduced which is 1 for the complete- 
form log-likelihood function and for the approximated one. 
The generalized expression for (1A.491 is shown in (1A.5H . For 
the complete-form log-likelihood function with 7 = 1, it is 
considerably hard to find a closed-form optimal solution to 
Qq : e due to the presence of [1 — cos(O gj e)]. Thus, 7 is set at 
which gives ( D~9) . To fill the gap between the complete and 
the approximated form of the log-likelihood function, a\ and 
a 2 H are replaced by |X 9) fc| 2 and \H q: k\ 2 respectively to take 
advantage of the modulation effect and frequency selectivity, 
which gives X 2 (q, k) in d2lT i. 



Appendix B 

Derivation of Variances and Biases in ([33]), OH), 

(|35), (OS) UNDER (A2) , (A3) , AND (A4) 

A. Bias Performance 

For DA, subtracting the left hand side of (|22) by the true 
value of <d q ,e yields 



A® q ,e = 6°A _ Q q e = arg Y, 



Xi(q,k)\ 2 (q,k) 



(B.52) 

Suppose that is in vicinity of O q ,e, tan(A0 9! e) can be 
approximated as 



tan(A0,, e ) = 



g [Ae,, fl ] 

3? [AQ 9i e] 



AG 



q .e 



(B.53) 



since AQ q g is sufficiently small. Expectation of AQ q g in 
(IB. 53b is approximated by 



Ae 



E{AM ~MwrMW} (B - 54) 

if E {5R [A9,, fl ]} > VVar {5R [A9 g , e ]}. Under (A2) , write 
E {3 [Ae qM ]} = -J2 sin(e g , e )A 2 (g, k) (B.55) 

E{5R[Ae g , e ]}= Y [cos(Q qt g) + SNR (q,k)]X 2 (q,k) 

(B.56) 



Under (A3) and (A4) , E {A6 9>0 } takes the form 
E {AG q , e } ~ Sin(6 - e) 



(B.57) 



cos(e 9 , e ) + SNR 

where SNR = SNR(<j,fc) and X 2 (q,k) = (2 + SNR- 1 ) = 
Const.. Assuming perfect phase unwrapping, ( |29l gives 
Bias{#DAj = Bias{V}. The 2x1 bias vector Bias{V} is 
associated with the Q x 1 bias vector Bias{V} by 



Bias{V} = (E 1 E)" 1 E J Bias{V} 



(B.58) 



where Bias{V} = 

3111(01,8) _ 

cos(ei, e )+SNR ' ' ' cosOn-i a)+SNR 

calculations lead to (OS) an d *0§)- 



B. Variance Performance 

Var {AO g e} can be formulated into 



sin(e 0- 8) 



<■ cos(e e)+SNR 
sinOq-x,,) _ 1T Standard 



Var { Ae, e} ^ Var{5[A9 ^l (B.59) 
(E{5R[Ae 9 , e ]}) 2 



under E {5R[A9 g , e ]} > ^Var {5i [AQ qy0 ]}. Again, assum- 
ing (A2) , (E {5R [A6 9 . e ]}) 2 can be written into (|B~60l . 

Note that, Var{3 [A6 ge ]} = E {(3 [AQ q e ]) 2 } - 
(E{3[A6 g , e ]}) 2 where 



(E{3[A8 g , e ]}) 2 = S in 2 (e g , e )^ J2 X ^ k ) 



(B.65) 
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e \ 2 )+23j{R<,,kR q ,N- k e- je i- e <t%<t%} 



f{Rq,klRq,N-k) 



( -^ 2 x^H+^w)(\^„,k\ 2 + \R q ,N~kl 

cxp V 2^^ CT ^(i- 7 co S (e g . a 7)+^;(i-7 2 ) 

^[20^.0^0^(1 - 7 cos(6 9 , e )) + - 7 2 )] 



(A.51) 



(\X q M 2 \H q ,k\ 2 + a 2 w cos(6 9 , e )) X 2 (q, k) 



(9[A9,, e ]) 2 = V V |^, fel | 2 |^, fe2 | 2 sm 2 (e 9)e )-|^ (? , fcl | 2 sin(e, ie )3{A3( g ,fc 2 )} 
— — _ V v y v .. ' 

kieZt fe 2 Glt 



^2 



|W g , fc2 | 2 sin(e 9 , e )9{A 3 (g,fc 1 )} + 3{A 3 ((Z,A ;i )}9{A3(g,fc 2 )} x A 2 (<z, fci)A 2 ( g , fe) 



Y 3 Y 4 Y 5 

X 3 (q,k) = e^- N ^X q ^H q ^W^ N . k +e^-^X qtN _ k H^ N _ k W q ^ 

E{ ]T E y i xF 5 } = sin 2 (6 g>e )o^ J] E A 2 (g,fci)A 2 (g,fc 2 ) 
fcj ei+ fe 2 gXj*" fei ezj" fc 2 exj 

E{ X] E x y s} « t^r (1 - cos(Q q ,e )) ^ \X q , k \ 2 \H q , k \ 2 \ 2 2 { qi k) 



(B.60) 



(B.61) 

(B.62) 
(B.63) 

(B.64) 



Expanding (3 [AOg.e]) 2 into the form shown as (IB.6U 
with X 3 (q,k) shown as (1B.621 1. E{£ fci£l + Sfe 2e x+ Y% x 
F5} and E {J2 kl ex+ £fc 2 ez+ ^4 x ¥5} can be computed 
as (lB~63l > and 'dfel, 'while E {£* iez + E fee i+ *a X 
Y 5 } and Eg^^^ x y s} can be neglected. 
Var {5 [A8 9 e]} under (A3) and (A4) is expressed by 

cos(6 9i e)) 



Var{3[Ae^]}«SJV SIXH " 



(2 + SNR- 1 ) 2 



(B.66) 



Substituting Var {3 [A6 9 . e ]} and (E {5i [A6, ifl ]}) 2 back into 
dE39l yields 



Var{A8 g , e } 



1 - cos(9 g ,e) 
SAT x SNR 



(B.67) 



Therefore, 

Var{^I} 



Var{V} = Var{(E T E)" 1 E T V} (B.68) 



Lengthy calculations lead to ( f33l > and d34l . Under (A3) and 
(A4) , the DA estimator reduces to NDA. 
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